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Abstract 

A theory for the effect of impurity scattering on electronic Raman scattering 
in unconventional superconductors is presented. The impurity dependence of 
the the spectra can be used to to distinguish between conventional (anisotropic 
s— wave) or unconventional (d-wave) energy gaps, and excellent agreement 
with the data on three cuprate superconductors is shown for (1^2 _y2 pairing. 
The full frequency dependence of the spectra can be reproduced by including 
inelastic scattering due to spin fluctuations and/or Coulomb scattering on a 
partially nested Fermi surface. 

PACS numbers: 78.30.-j, 74.72.-h, 74.25. Gz, 74.20.Mn 



1 



Typeset using REVT^ 



By now a more coherent picture of the pair state symmetry in cuprate superconductors 
has emerged, and it has become increasingly evident that a gap with predominantly d^2_y2 
pairing symmetry can provide an adequate description of many experiments However, 
the question remains or not whether the order parameter changes sign around the Fermi 
surface. The available data on Josephson junction and point-contact tunneling measure- 
ments seems to in favor of the former, although examples exist to the contrary In 
addition, it has been argued that the impurity dependence of certain correlation functions 
can be used to determine whether the gap changes sign. In particular, a detailed analysis 
has been given for the case of the electromagnetic penetration depth |^, and the crossover 
from T to at a temperature T* which grows with impurity concentration has be taken 
to support (i— wave pairing. Similar considerations can be made for the Knight shift, NMR 
rates, etc. However, the evidence gained from the penetration depth or other correlation 
functions is incomplete in that only the topology of the energy gap nodes can be identified 
from the power-law behavior at low temperatures, which leaves the actual gap symmetry 
unspecified. For instance, in tetragonal superconductors, there are 5 pure representations, 
{dx2-y2,drcy,d3r2-i,dxz sud dyz), which yield the same low temperature power-laws and im- 
purity dependences for, e.g., the penetration depth, NMR rates, density of states, etc. 

It has recently been shown that both the electronic and phonon contributions to Raman 
scattering in superconductors can be used to identify the order parameter symmetry and can 
distinguish between the various representations for d— wave pairing [^]. Raman scattering 
measures effective intracell density fluctuations Pq, 

Pq = E 7(k)ci_q/2,,Ck+q/2,<T, (l) 

No- 
where cj^o.(ck,(T) is the creation (annihilation) operator for an electron with wavevector k and 
spin a. For the case of non-resonant scattering 0, the Raman vertex 7 is directly related 
to the curvature of the band dispersion e(k) and incident and scattered polarization 
light vectors via 
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which can be expanded in terms of a complete set of orthogonal functions on the Fermi 
surface, 7(k) = Y^lIl^lO^)- For bands with non-parabolic dispersion, all L symmetry 
channels can in principle contribute to the vertex. Thereby, fluctuations on different parts 
of the Fermi surface can be selected by orienting the polarization light vectors. 

It was shown in Ref. |^ that the symmetry- (polarization-) dependence of the Raman 
spectra straightforwardly arises through the coupling of the Raman vertex to an anisotropic 
energy gap. The main features of the theory for clean superconductors with d^'^-y'^ pairing 
symmetry are as follows: 1) due to the presence of the nodes, no clear gap threshold is 
seen for any polarization channel, 2) the peak in the spectra occurs at the largest frequency 
for polarizations which select the Big channel (i.e, for YBCO, crossed in-lane polarizations 
aligned 45 degrees with respect to the a — h axes), while lower peaks are observed in the 
(45 degree rotation of the Big orientation) and Aig (mostly parallel polarizations) channels 
in descending order, respectively, and 3) the low frequency behavior is given by uj'^ for the Big 
channel, while is linear in uj for the other channels. Such a strong polarization dependence 
of the spectra allowed for a stringent test of the predictions for dx2_y2 pairing, and good 
agreement of the theory was shown with the data on three cuprate materials, YBa2Cu307, 
Bi2Sr2CaCu208, and Tl2Ba2Cu06. 

However, like other correlation functions, only | A(k) | could be determined from the 
fit, and thus it could not be ascertained whether the gap changed sign around the Fermi 
surface. In addition, the fit was less satisfactory for samples with increasing disorder, with 
the best agreement found for a sample of Bi 2:2:1:2 with a transition temperature ~ 86 K. 
The theory also failed to predict any Raman intensity in the normal state due to phase space 
restrictions (consequence of the limit q ^ 0), and a large and unspecified smearing width 
F ~ 0.2Ao was needed even at low temperatures. 

The purpose of this Letter is to show that impurity scattering can resolve the above 
deficiencies of the theory. The T— matrix approach to incorporate repeated scattering events 
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from a single impurity site has been well documented 0. The two parameters to describe 
the theory are the cotangent of the scattering phase shift, c = cot(5), and the impurity 
concentration rii described through the scattering rate F = rii/nNp, where Np is the density 
of states per spin at the Fermi level. To maintain gauge invariance, the effects of impurity 
scattering are included into both the self energy and vertex renormalizations. The case 
for an even parity vertex (such as the Raman vertex) has been worked out for general pair- 
states under the assumption of particle-hole symmetry and s-wave impurity scattering in Ref. 
P], and thus we can adopt these methods accordingly. For the case of the Raman vertex, 
impurity dressings are only important for the case of the density (L = 0) channel, while other 
channels do not couple to the scattering. Since the density channels are completely screened 
out (intercell charge fluctuations) for vanishing momentum transfers, for our purposes we can 
set the vertex corrections to zero. This will not hold for an impurity interaction which is non- 
zero in other channels (i.e., a k-dependent impurity interaction). This can be remedied by 
including couplings to other channels, but the analysis can not be carried as far analytically. 
Our expression for the Raman response is then given by 

iuj k 

where Tr denotes the trace, f are Pauli matrices, and 7 is the bare Raman vertex, Eq. 
(2). We use the standard weak coupling BCS Green's function dressed by the impurity self 
energy, 

^^'^'^^)=(z^)2-e(k)^-A(k)2' 
where the tilde indicates the renormalized frequency and band energy via 

iCb = iuo — So(itD), 
e(k) = e(k)-S3(z(:;). (5) 

Lastly, the matrix self energy is given in terms of the integrated Green's function g{iuj) = 
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Our choice of the gap parameter is A(k) = Aq cos (20), where we will for simplicity chose to 
work with a 2-D (cylindrical) Fermi surface. 

The results of the theory with the inclusion of impurity scattering in the Born (c >> 1) 
and unitary (c = 0) limits are shown in Figure 1. The impurity scattering smears out the 
sharp features of the spectra (i.e., the logarithmic divergence of the spectra at the gap edge 
for the Big channel) and thus can provide an explanation for the broadening seen in the data 
at low temperatures. In addition, the anisotropy of the peak positions is maintained (unless 
the scattering is so large that the normal state behavior is recovered). The coupling of the gap 
and Raman vertex lead to polarization dependent peak positions, with Upeak/^o = 2.0, 1.8, 
and 1.2 for the Big,B2g, and Aig channels, respectively. Thereby, the channel which shows 
the peak at the highest frequency gives the predominant symmetry of the energy gap. This 
was used in Ref. ^ to conclude that the energy gap in the cuprates is predominantly of 
dx2-y2 symmetry as opposed to dxy or other c?— wave representations. 

To determine further whether the gap is entirely of dx2_y2 symmetry, the impurity depen- 
dence of the low frequency behavior of the channel dependent spectra can be exploited. From 
Figure 1, it is seen that Born impurity scattering yields low frequency exponents that are 
the same as for clean materials (unless of course F ~ Aq and the normal state is recovered). 
However, in the unitary case, while the low frequency behavior remains linear in frequency 
in the Aig and B2g channels, below a characteristic frequency uo* the behavior changes from 
uo^ to linear in uj for the Big channel. This is due to a nonzero density of states at the Fermi 
level, which allows for normal-state-like behavior to be recovered 0. As in the case of the 
penetration depth , the scale uj* grows with increasing impurity concentrations. However, 
the exponent is symmetry dependent (remains 1 for Aig and i?2g channels, while decreases 
from 3 to 1 for the Big channel, which is unlike the penetration depth). This is in marked 
contrast to the impurity dependence of the spectra for anisotropic s-wave gaps, since in 
that case, the exponents for all channels would grow as the gap becomes more isotropic for 



increased impurity scattering . In particular, the impurity dependence of the Big channel 
is opposite to what one would expect if the gap was anisotropic s-wave. Thus the impurity 
dependence can be systematically checked to determine whether the gap has accidental or 
intrinsic zeroes, or more generally, if the gap is anisotropic s-wave with predominantly Big 
symmetry, or if it is dx2_y2 symmetry. 

To make a comparison to experiment, we note that the impurity dependence, for in- 
stance via Zn doping, has not been systematically checked via Raman scattering for any of 
the cuprate superconductors. Raman scattering measurements of the electronic background 
have been performed on materials of varying quality. In particular, we focus on three mate- 
rials, as-grown Bi2Sr2CaCu208 (Tc = 86/^), oxygen-annealed Bi2Sr2CaCu208 (Tc = 
0, and Tl2Ba2Cu06 (Tc = SOfC) [1^. The materials show polarization dependent spectra, 
with the peak positions accurately predicted using a gap of dx-^_y2 symmetry (see Ref. 0]). 
These compounds differ with respect to both surface quality and the presence of bulk im- 
purity defects. Fits using the theory for clean materials |^ required a Gaussian broadening 
width r ~ 0.2Ao. Such large broadening, if attributed solely to impurity scattering, would 
of course also lead to a large reduction of T^. However, the majority of the smearing could 
be due to a distribution of Tc values (of the order of 20 percent in Tl 2:2:0:1) and inelastic 
quasiparticle scattering, which could be taken into account via strong coupling theory. The 
actual magnitude of the scattering to be attributed to impurities remains an open issue 
which requires both a non-phenomenological theory to include the effects of inelastic scat- 
tering, doping, and disorder on the pairing interactions, as well as a systematic experimental 
analysis. 

We present the comparison of the theory to the spectra obtained in Refs. P,p!0| on the 



three compounds in Figure 2. The spectra shown in the Figure are for the Big and i?2g 
channels for both Bi 2:2:1:2 samples, and for Big and Aig for Tl 2:2:0:1. The parameters 
used to fit the data are Aq = 287, 192, and 240 cm~^, for the energy gap, and resonant 
impurity scattering rate F/Aq = 0.125, 0.2, and 0.25 for as-grown, oxygen- annealed Bi 
2:2:1:2, and Tl 2:2:0:1, respectively. No Gaussian broadening was used. The theory shows 



that the anisotropy of the peak positions is preserved even in the presence of impurity 
scattering which suggests that the gap anisotropy is robust to disorder. The hnear rise 
in frequency of both the Big and Aig channels is a direct consequence of resonant impurity 
scattering. This is shown in log-log plot of the low frequency part of the Big response for the 
three materials in the inset of Figure 2. Born scattering can not fit the data as the amount of 
scattering to produce the linear rise in the Big channel fully suppresses the peak feature. The 
crossover scale u* grows with increasing disorder, which leads to a more pronounced linear 
rise of the Big channel. The crossover frequency is given by u* /Aq = 0.38, 0.47, and 0.73 for 
the as-grown, oxygen-annealed Bi 2:2:1:2, and Tl 2:2:0:1, respectively. This is in consistent 
with the predictions of rf— wave pairing and is difficult to reconcile using an anisotropic 
s— wave gap. Thus, this crucial issue merits a systematic experimental examination of the 
low frequency region of the spectra using greater resolution. 

The theoretically predicted spectra fail to reproduce the flat background extending out 
to large Raman shifts [|T2| since the calculated spectra falls off as l/uj [O]. "Marginal-like" 



behavior can be phenomenologically modelled by inserting a self energy piece of the 
form 



E"{u;) = a^u;^ + PT^, (7) 

where a and /? are constants. This form for the self-energy arises in the normal state when 
electron-electron scattering on a nested Fermi surface is taken into account [0 or spin- 



fluctuation scattering is considered [|T6|. While the inclusion of this term to the self energy 
does not take the modifications of superconductivity into account and also neglects vertex 
corrections of the interactions, the flat behavior at high frequencies can be reproduced to 
flt the data. This is shown explicitly in Figure 3, which compares the results of the theory 
with and without the added self energy to the Big data on the annealed Bi 2:2:1:2 sample 
using a value of the gap, Aq = 190 cm^^. The parameters used are a/Ao = 0.25, /5 = 3.3, 
and small resonant scattering F = O.OIAq for the solid line, which is to be compared to the 
flt made in Figure 2 without the self energy term added (redrawn in Figure 3 as the dotted 



line). While the low frequency portion of the spectrum uj < 2Ao is relatively unaffected, 
the higher frequency region is well fit by the self energy inclusion. It is expected that a 
better fit could be obtained when modifications due to superconductivity are included in 
the interaction self energy. 

The agreement of the theory to the data leads to the following conclusions: First, the 
anisotropy of the peak positions result from a coupling of the Raman vertex to an energy gap 
with predominantly d^-i-y^ symmetry. Second, since the spectra do not show evidence for a 
reorganization of states to higher frequencies with increasing disorder, and on the contrary, 
reorganize to lower frequencies, this strongly suggests that the gap therefore changes sign 
and has unconventional dj.2_y2 symmetry as opposed to conventional anisotropic s-wave 
symmetry. This is seen from the increased prominence of the linear rise of the Big spectra 
and the increase of uj* for increasing resonant scattering. Lastly, the full frequency range 
of the spectra can be modelled by taking into account an additional interaction self energy. 
A systematic check of impurity scattering on a single compound would clarify these points 
further. 
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FIGURES 

FIG. 1. Raman spectra for Big,B2g, and Aig channels, with values of c and F indicated. 
Magnitude of the Raman vertices are set to 1. The inset shows the low frequency behavior of the 
Big channel. For the case of resonant scatterers (c = 0, bottom panel), the low frequency behavior 
crosses over from w to w^. 

FIG. 2. Fit to the data taken for the Big and B2g channels in as- grown Bi 2:2:1:2 (Tc = 86K, 
top panel), O2 annealed Bi 2:2:1:2 (Tc = 79K, middle panel) and for the Big and Aig channels in 
Tl 2:2:0:1 (Tc = 80K, bottom panel) obtained in Refs. and |jl^. The phonon contributions have 
been subtracted out in the middle and bottom panels. The parameters used are Aq = 287, 190, and 
240 cm~^, F/Aq = 0.125, 0.2, and 0.25, for the top, middle, and bottom panels, respectively. 
Inset: Log-log plot of the low frequency portion of the Big response. The crossover frequency 
uj* /Aq = 0.38,0.47, and 0.73 for the top, middle, and bottom panels, respectively. 

FIG. 3. Theory curves calculated with and without interaction self energy, Eq. 7, compared to 
the Raman spectrum for the Big channel of the O2- annealed Bi 2:2:1:2 sample. The dotted line 
is redrawn from Fig. 2. 
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